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ABSTRACT 

We report the discovery of four relatively massive (2-7 Mj) transiting extrasolar planets. HAT-P- 
20b orbits the moderately bright V=11.339 K3 dwarf star GSC 1910-00239 on a circular orbit, with a 
period P = 2.875317 ±0.000004d, transit epoch = 2455080.92661 ±0.00021 (BJDutc), and transit 
duration 0.0770±0.0008d. The host star has a mass of 0.76±0.03Mo, radius of 0.69±0.02 Rq, effective 
temperature 4595±80K, and metallicity [Fe/H] ~ +0.35±0.08. The planetary companion has a mass 
of 7.246±0.187A/j, and radius of 0.867±0.033iT!j yielding a mean density of 13.78±1.50g cm'^, which 
is the second highest value among all known exoplanets. HAT-P-21b orbits the V=11.685 G3 dwarf 
star GSC 3013-01229 on an eccentric (e = 0.228±0.016) orbit, with a period P = 4.124481±0.000007d, 
transit epoch = 2454996.41312 ±0.00069, and transit duration 0.1530 ±0.0027d. The host star has 
a mass of 0.95 ± 0.04 M©, radius of 1.10 ± 0.08 Rq, effective temperature 5588 ± 80 K, and metallicity 
[Fe/H] = +0.01 ± 0.08. The planetary companion has a mass of 4.063 ± 0.161 Mj, and radius of 
1.024 ± 0.092 i?j yielding a mean density of 4.68to Jg g cm'^. HAT-P-21b is a border-line object 
between the pM and pL class planets, and the transits occur near apastron. HAT-P-22b orbits the 
bright V=9.732 G5 dwarf star HD 233731 on a circular orbit, with a period P = 3.212220±0.000009d, 
transit epoch = 2454930.22001 ± 0.00025, and transit duration 0.1196 ± 0.0014 d. The host star 
has a mass of 0.92 ± 0.03 M©, radius of 1.04 ± 0.04 7?©, effective temperature 5302 ± 80 K, and 
metallicity [Fe/H] = +0.24 ± 0.08. The planet has a mass of 2.147 ± 0.061 Afj, and compact radius 
of 1.080 ± 0.058 i?j yielding a mean density of 2.1ltg j° gem ^. The host star also harbors an M- 
dwarf companion at a wide separation. Finally, HAT-P-23b orbits the V=12.432 GO dwarf star 
GSC 1632-01396 on a close to circular orbit, with a period P ^ 1.212884 ± 0.000002 d, transit epoch 
Tc = 2454852.26464 ± 0.00018, and transit duration 0.0908 ± 0.0007 d. The host star has a mass of 
1.13±0.04M©, radius of 1.20 ± 0.07 i?©, effective temperature 5905±80K, and metallicity [Fe/H] = 
+0.15 + 0.04. The planetary companion has a mass of 2.090 + 0.111 Mj, and radius of 1.368 + 0.090 i?j 
yielding a mean density of 1.01 + 0.18 gcm~'^. HAT-P-23b is an inflated and massive hot Jupiter on 
a very short period orbit, and has one of the shortest characteristic in-fall times (7.5l^ g Myr) before 
it gets engulfed by the star. 

Subject headings: planetary systems ~ stars: individual ( HAT-P-20, GSC 1910-00239, HAT-P-21, 
GSC 3013-01229, HAT-P-22, HD 233731, HAT-P-23, GSC 1632-01396 ) techniques: 
spectroscopic, photometric 
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The majority of the ~ 90 known transiting extraso- 
lar planets (TEPs) have been found to lie in the 0.5 Mj 
to 2.0 Mj mass range. The apparent drop in their 
mass distribution . at_~2Mj ha s been noted by, e.g., 
iSouthworth et al.l (|2009D . and bv lTorres et"al] (|2010 ^. In 
the currently known sample, 75% of the TEPs have plan- 
etary mass Mp < 2.0 Mj, and there appears to be a mi- 
nor peak in their occurrence rate at Mp « 2Mj, which 
then sharply falls off towards higher masses. Are there 
any biases present against discovering massive plan- 
ets? Such planets tend to be less inflated, and the- 
ory dictates that their radii shrink as their mass in- 
cr eases towards the brown dwarf reg ime. According 
to iBarafFe. Chabrier. fc BarmanI (|2010( 1 , this reversal of 
the Mp — Rp relation happens around Mp « 2 — 3Mj, 

and falls off as Rp oc Mp^^^ (see e.g. iFortnev et al.l 
l2009f ). The smaller radii for massive planets yield a 
minor bias against discovering them via the transit 
method, since they produce shallower transits. Very mas- 
sive p lanets can induce st ellar variability of their host 
stars (|Shkolnik et al.|[2009() , somewhat decreasing the ef- 
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ficiency of detecting their shallow transits via simple al- 
gorithms that expect constant out-of-transit light curves. 
Also, the host stars of massive planets are typically more 
rapid rotators: the average vsini for host stars with 
planets Mp < 2Mj is 3.9 ms~^ (with 3.8 ms~^ standard 
deviation around the mean), whereas the same values 
for the massive planet hosts stars, are 10.9ms~^(with 
12.7 ms~^ standard deviation around the meanl^. The 
five fastest rotators all harbor planets more massive than 
2Afj. This presents a bias against discovering them ei- 
ther via radial velocity (RV) searches, which are more 
efficient around quiet non-rotating dwarfs, or via transit 
searches, where the targets may be discarded during the 
confirmation phase. Along the same lines, the large RV 
amplitude of the host star, as caused by the planetary 
companion, may even lead to erroneous rejection dur- 
ing the reconnaissance phase of candidate confirmation, 
since such systems resemble eclipsing binaries. Finally, 
there is a tenden cy that massive planets are more likely 
to be eccentricI3 (jSouthworth et al]|2009D , meaning that 
they require more RV observations for proper mapping of 
their orbits, and thus leading to a slower announcement 
rate. On the other hand, a strong bias for detecting such 
planets-compensating for most of the effects above-is the 
fact that the large RV amplitudes of the host stars are 
easier to detect, since they do not require internal pre- 
cisions at the ms"^ level (see HAT-P-2, where valuable 
data was contributed to the RV fit by mod est precision 
instruments yielding ~ lkms~^ precision; iBakos et al.l 
|2007| ). Altogether, while there are minor biases for and 
against detecting massive transiting planets, their overall 
effect appears to be negligible, and the drop in frequency 
at > 2 Mj seems to be real. 

Massive planets are important for many reasons. They 
provide very strong constraints on formation and migra- 
tion theories, which need to explain the observed distri- 
bution of p lanetary system parameters in a wide range 
(|Baraffe et al. 2008:i|Barafer Chabricr, & Barman 2010), 
from 0.01 Afj (Corot-7b; Quclo z et al.ll2009l) to 26.4 Afj 
(Corot-3b; Deleuil et al. 2008,). Heavy mass objects 
necessitate the inclusion of other physical mechanisms 
for the for mation and migration, such as planet-plane t 



scattering ()Chatteriee et all l2008t iFord fc Rasiol 



and the Kozai-mechanism (iFa brvckv &: Tremaind 
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They are border-line objects between planets and brown- 
dwarfs, and help us u nderstand how thes e populations 
differ and overlap (see iLeconte et al.l[2009l for a review) . 
For example, a traditional definition of planets is that 
they have no Deuterium burning, where the Deuterium 
burning limit is thought to be around 13Afj. However, 
there are large uncertainties on this limit due to the nu- 
merous model parameters and solutions, and the fact 
that Deuterium may be able to burn in the H/He lay- 
ers above the core (|Baraffe et al.l[2008[ ) . Another possible 
definition of planets is based on their formation scenario, 
i.e. they are formed by accretion in a protoplanetary disk 
around their young host star, as opposed to the gravita- 
tional collapse of a molecular cloud (brown dwarfs). 

Perhaps related to the formation and migration 
mechanisms, a number of interesting correlations 
involving massive planets have been pointed out. 

This includes those 4 planets announced in this paper. 
See e.g. exoplanets.org for statistics. 



lUdrv et all ()2002D noted that short period massive plan- 
ets are predominantly f ound in binary stellar systems. 
iSouthworth et al.l ()2009[ ) noted that only 8.6% of the low 
mass planets show significantly eccentric orbits, whereas 
77% of the massive planets have eccentric orbits (al- 
though low- mass systems have lower S/N RV curves, ren- 
dering the detection of eccentric orbits more difficult). 
Curiously, there appears to be a lack of correlation be- 
tween planetary mass and host star metallicity, while one 
would naively think that the formation of high mass plan- 
ets (via core accretion) would require higher metal con- 
tent. Until this work, there was a hint of a correlation 
between planetary and stellar mass (e.g. iDeleuil et al.l 
[2008I) . in the sense that the most massive planets orbited 
> 1.2 M0 stars, and there was a (biased) tendency 
that lower mass planets orbit less massive stars. 

All of these observations suffer from small-number 
statistics and heavy biases. One way of improving our 
knowledge is to expand the sample of well-characterized 
planets. In this work we report on 4 new massive tran- 
siting planets around bright stars, namely GSC 1910- 
00239, GSC 3013-01229, HD 233731, and GSC 1632- 
01396. This extends the currently known sample of 
bright {V < 13.5) and massive {Mp > 2 Mj) transiting 
planets by 30% (from 13 to 17). These discoveries were 
made by the Hu ngarian-made Auto mated Telescope Net- 
work (HATNet; IBakos et al.ll2004[ ) survey. HATNet has 
been one of the main contributors to the discovery of 
TEPs, among others s uch as the g round-based Super- 
WASP (IPollacco et al.l[2006iV TrES (lAlonso et al.l [200l 
and XO projects ( McCullough et al. 20051) . and s pace- 
borne searches such as CoRoT m etall [20061 ) and 
Kepler (jBorucki et al.l [20Tol) . In operation since 2003, 
HATNet has now covered approximately 14% of the sky, 
searching for TEPs around bright stars (8 < / < 14). 
We operate six wide-field instruments: four at the Fred 
Lawrence Whipple Observatory (FLWO) in Arizona, and 
two on the roof of the hangar servicing the Smithso- 
nian Astrophysical Observatory's Submillimeter Array, 
in Hawaii. 

The layout of the paper is as follows. In Section [2] we 
report the detections of the photometric signals and the 
follow-up spectroscopic and photometric observations for 
each of the planets. In Section [3] we describe the anal- 
ysis of the data, beginning with the determination of 
the stellar parameters, continuing with a discussion of 
the methods used to rule out nonplanetary, false posi- 
tive scenarios which could mimic the photometric and 
spectroscopic observations, and finishing with a descrip- 
tion of our global modeling of the photometry and radial 
velocities. Our findings are discussed in Section |4l 

2. OBSERVATIONS 

2.1. Photometric detection 

Table [T] summarizes the HATNet discovery observa- 
tions of each new planetary system. The calibration 
of the HATNet frames was carried out using standard 
procedures correcting for the CCD bias, dark-current 
and fiatfield structure. The calibrated images were 
then subje cted to star d e tectio n and astrometry, as de- 
scribed in iPal &: BakosI (|2OO60 . Aperture photometry 
was performed on each image at the stellar centroids 
derived from the Two Micron All Sky Survey (2MASS; 
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iSkrutskie et alj [20061 ) catalog and the individual astro- 
metric solutions. For certain datasets (HAT-P-20, HAT- 
P-22, H AT-P-23) we also carried out an image sub- 
traction (|Alardl [20000 based photo metric reduction us- 
ing disc rete kernels ( Bramich[ 120081 ). as described in iPall 
(|2009b( ). The resulting light curves were decorrelated 
(cleaned of trends) using the Externa l Parameter Decor- 
relation (EPD; see lBakos et al.|[2010[ ) technique in "con- 
sta nt" mode and the Trend Filtering Algorithm (TFA; 
see lKovacs et al.l [20051 ) . The light curves were searched 
for periodic box-s haped signals usin g the Box Least- 
Squares (BLS; see lKovacs et al.l '2002) method. We de- 
tected significant signals in the light curves of the stars 
as summarized below: 

• HAT-P-20 - GSC 1910-00239 (also known as 
2MASS 07273995-f2420118; a = 07h27"^39.96", 
6 = -f24°20'11.9"; J2000; V=11.339). A signal 
was detected for this star with an apparent depth 
of ~10.8mmag, and a period of P =2.8753 days 
(see Figure [Ij . Note that the depth was attenuated 
by the presence of the fainter neighbor star that is 
not resolved on the coarse resolution (14"pixel~^) 
HATNet pixels. Also, the d epth by fitting 
a trap ese instead of the correct iMandel fc Agoll 
()2002f ) model, is so mewhat shallower than the max- 
imum depth in the lMandel &: Agoll (|2002D model fit 
(which was 19.6 mmag; see later in § 13. 3p . The drop 
in brightness had a first-to-last-contact duration, 
relative to the total period, of q =0.0268 ± 0.0003, 
corresponding to a total duration of Pq =1.848 ± 
0.019 hr. HAT-P-20 has a red companion (2MASS 
07273963-1-2420171, J - K = 0.92) at 6.86" sep- 
aration that is fainter than HAT-P-20 by Ai? = 
1.36 mag. 

• HAT-P-21 - GSC 3013-01229 (also known as 
2MASS 11250598-f4101406; a = llh25'"05.88", 
5 = -f41°01'40.6"; J2000; V=11.685). A signal 
was detected for this star with an apparent depth 
of '^8.0mmag, and a period of P =4.1245 days (see 
Figure [2|) . The drop in brightness had a first-to- 
last-contact duration, relative to the total period, 
of q =0.0371 ±0.0007, corresponding to a total du- 
ration of Pq =3.672 ± 0.065 hr. 

• HAT-P-22 - HD 233731 (also known as 
GSC 03441-00925 and 2MASS 10224361+5007420; 
a = 10h22™43.68^ 5 = -H50°07'42.0"; J2000; 
V=9.732). A signal was detected for this star with 
an apparent depth of ~9.7mmag, and a period of 
P =3.2122 days (see Figure [3]). The drop in bright- 
ness had a first-to-last-contact duration, relative to 
the total period, of q =0.0372±0.0004, correspond- 
ing to a total duration of Pq =2.869 ± 0.033 hr. 
HAT-P-22 has a close red companion star (2MASS 
10224397+5007504, J - K = 0.86) at 9.1" 
separation and Ai = 2.58 magnitude fainter. 

• HAT-P-23 - GSC 1632-01396 (also known as 
2MASS 20242972+1645437; a = 20h24'"29.88", 
S = +16°45'43.7"; J2000; V=12.432). A signal was 
detected for this star with an apparent depth of 
~ 11.5 mmag, and a period of P =1.2129 days (see 



Figure [J). Similarly to HAT-P-20, the depth was 
attenuated by close-by faint neighbors. The drop 
in brightness had a first-to-last-contact duration, 
relative to the total period, of q =0.0748 ± 0.0006, 
corresponding to a total duration of Pq =2.178 ± 
0.017 hr. 
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Fig. 1.— Unbinned light curve of HAT-P-20 including all 2600 
instrumental R band 5.5 minute cadence measurements obtained 
with the HAT-7 and HAT-8 telescopes of HATNet (see Table[T]for 
details), and folded with the period P = 2.8753172 days resulting 
from the global fit described in Section [3] The solid line shows the 
"PIPS" transit model fit to the light curve (Section l3.3|l . The bold 
points in the lower panel show the light curve binned in phase with 
a bin-size of 0.002. 
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Fig. 2.— Unbinned hght curve of HAT-P-21 including all 28,000 
instrumental I band and R band 5.5 minute cadence measurements 
obtained with the HAT-5, HAT-6, HAT-8 and HAT-9 telescopes of 
HATNet (see Tabled] for details), and folded with the period P = 
4.1244808 days resulting from the global fit described in Section |3] 
The soli d lin e shows the "PIPS" transit model fit to the light curve 
(Section lS.3| l. The bold points in the lower panel show the light 
curve binned in phase with a bin-size of 0.002. 



2.2. Reconnaissance Spectroscopy 

As is routine in the HATNet project, all candidates 
are subjected to careful scrutiny before investing valu- 
able time on large telescopes. This includes spectro- 
scopic observations at relatively modest facilities to es- 
tablish whether the transit-like feature in the light curve 
of a candidate might be due to astrophysical phenom- 
ena other than a planet transiting a star. Many of these 
false positives are associated with large radial-velocity 
variations in the star (tens of kms~^) that are easily rec- 
ognized. The reconnaissance spectroscopic observations 
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TABLE 1 

Summary of photometric observations 



Instrument /Field 
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^ Observations were performed without guiding due to a technical problem with the guiding system, and 
resulted in decreased data quality. 
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Fig. 3.— Unbinned light curve of HAT-P-22 including all 4200 
instrumental R band 5.5 minute cadence measurements obtained 
with the HAT-5 telescope of HATNet (see the text for details), 
and folded with the period P = 3.2122198 days resulting from the 
global fit described in Section [3] The solid l ine shows the "P1P3" 
transit model fit to the light curve (Section 13.31 1. The bold points 
in the lower panel show the light curve binned in phase with a 
bin-size of 0.002. 

and results for each system are summarized in Table 121 
below we provide a brief description of the instruments 
used, the data reduction, and the analysis procedure. 

One of the tools we have used for this purpose is the 
Harvard-Smithsonian C enter for Astro physics (CfA) Dig- 
ital Speedometer (DS; iLathaml Il992j ). an echelle spec- 
trograph mounted on the FLWO 1.5 m telescope. This 
instrument delivers high- resolution spectra (A/ A A « 
35,000) over a single order centered on the Mg lb triplet 
(~5187A), with typically low signal-to-noise (S/N) ra- 
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Fig. 4.— Unbinned light curve of HAT-P-23 including all 3500 
instrumental R band 5.5 minute cadence measurements obtained 
with the HAT-6 and HAT-9 telescopes of HATNet (see Table [T] for 
details), and folded with the period P = 1.2128841 days resulting 
from the global fit described in Section \3\ The soli d line shows the 
"P1P3" transit model fit to the light curve (Section 13.31 1. The bold 
points in the lower panel show the light curve binned in phase with 
a bin-size of 0.002. 

tios that are nevertheless sufficient to derive radial veloc- 
ities (RVs) with moderate precisions of 0.5-1.0kms~^ for 
slowly rotating stars. The same spectra can be used to 
estimate the effective temperature, surface gravity, and 
pr ojected rotati o nal ve locity of the host star, as described 
bv lTorres et al.l (|2002| ). With this facility we are able to 
reject many types of false positives, such as F dwarfs or- 
bited by M dwarfs, grazing eclipsing binaries, or triple or 
quadruple star systems. Additional tests are performed 
with other spectroscopic material described in the next 
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section. 

Another of the tools we have used for this purpose is 
the Fibre-fed Echelle Spectrograph (FIES) at the 2.5 m 
Nordic Optical Telescope ( NOT) at La Palma, Spain 
([Diupvik &: AndersenI I2010D . We used the medium- 
resolution fiber which produces spectra at a resolution 
of A/ A A « 46,000 and a wavelength coverage of ~ 3600- 
7400 A to observe HAT-P-20. The spectrum was ex- 
tracted and analyzed to measure the radial velocity, ef- 
fective temperature, surface gravity, and projected rota- 
tion velocity of the host star, follo wing the procedures 
described by iBuchhave et al.l ()2010[) . 

Based on the observations summarized in Table [2] 
we find that HAT-P-21, HAT-P-22 and HAT-P-23 have 
rms residuals consistent with no detectable RV variation 
within the precision of the measurements. Curiously, 
HAT-P-20 showed significant RV variations, even at the 
modest (~ 0.5kms~^) precision of the Digital Speedome- 
ter, and the reconnaissance RV variations (including the 
FIES spectrum; see later) phased up with the photomet- 
ric ephemeris. All spectra were single- lined, i.e., there 
is no evidence that any of these targets consist of more 
than one star. The gravities for all of the stars indicate 
that they are dwarfs. 

2.3. High resolution, high S/N spectroscopy 

We proceeded with the follow-up of each candidate by 
obtaining high-resolution, high-S/N spectra to character- 
ize the RV variations, and to refine the determination of 
the stellar parameters. These observations are summa- 
rized in Table [31 The RV measurements and uncertain- 
ties are given in Table S Tabled Table Hand Table [7] for 
HAT-P-20 through HAT-P-23, respectively. The period- 
folded data, along with our best fits described below in 
Section [31 are displayed in Figure [Bl through Figure [SI for 
HAT-P-20 through HAT-P-23. Below we briefiy describe 
the instruments used, the data reduction, and the anal- 
ysis procedure. 




Fig. 5. — Keck/HIRES guider camera snapshots of HAT-P-20 
through HAT-P-23 (labeled). North is up and East is to the left. 
The snapshots cover an area of approximately 30 X 20". The slit 
is also visible, as positioned on the planet host stars. 

Observations were made of all four planet h ost stars 
with the HIRES instrument ()Vogt et al.l I1994D on the 
Keck I telescope located on Mauna Kea, Hawaii. The 
width of the spectrometer slit was 0'.'86, resulting in a re- 
solving power of A/AA « 55,000, with a wavelength cov- 



erage of -3800-8000 A. We typically used the B5 decker 
yielding a 3.5"{H) x 0.861"(M^) slit, and for the last few 
observations on each target we used the C2 decker that 
enables a better sky subtraction due to the longer slit 
14.0"(iJ) X 0.861"(W^). The slit height was oriented with 
altitude (vertical), except for rare cases, when the slit 
would have run through the faint companion to HAT- 
P-20 or HAT-P-22. A Keck/HIRES snapshot for each 
planet host star is shown in Figure [5l Spectra were ob- 
tained through an iodine gas absorption cell, which was 
used to superimpose a dense forest of I2 lines on the 
stellar spect rum and establish an accurate wavelength 
fiducial (see iMarcv fc Butleii I1992D . For each target an 
additional exposure was taken without the iodine cell, for 
use as a template in the reductions. Relative RVs in the 
sol ar system bar y centri c frame were derived as described 
by iButler et al.l 1)199 6') . incorporating full modeling of 
the spatial and temporal variations of the instrumental 
profile. 
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Fig. 6. — Top panel: Keck/HIRES RV measurements for 
HAT-P-20 shown as a function o f or bital phase, along with our 
best-fit eccentric model (see Table [T4t . Zero phase corresponds to 
the time of mid-transit. The center-of-mass velocity has been sub- 
tracted. The rms around the best orbital fit is 16.2ms~^. Second 
panel: Velocity O — C residuals from the best fit. The error bars 
for both the top and second panel include a component from the 
jitter ( IS.T ms"-*-) added in quadrature to the formal errors (see 
Section l3.3l l. Third panel: Bisector spans (BS), with the mean 
value subtracted. The me asurement from the template spectrum is 
included (see Section |3.2|I . Bottom panel: Relative chromospheric 
activity index S measured from the Keck spectra. 

In each of Figures [nHHl we show also the relative S in- 
dex, which is a measure of the chromospheric activity of 
the star derived from the fiux in the cores of the Ca II H 
and K lines. This index was computed following the pre- 
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TABLE 2 

Summary of reconnaissance spectroscopy observations 



Instrument 


Dato(s) 


Number of Spectra 


(K) 


log g* 
(cgs) 


V sin j 
(kms 


7rv'' 
(kms '-) 


HAT-P-20 














DS 
FIES 


2009 Feb 11-2009 Feb 15 
2009 Oct 07 


3 
1 


4500 ± 125 
4500 


4.0 ±0.25 
4.0 


0±4 
4 


-18.81 ±0.68 
-16.76 ±0.1 


HAT-P-21 














DS 


2009 Mar 08-2009 Apr 05 


3 


5750 ± 125 


4.5 ±0.25 


3.0 ±3.0 


-53.19 ±0.09 


HAT-P-22 














DS 


2009 Feb 11-2009 Feb 16 


4 


5250 ± 125 


4.5 ±0.25 


2.0 ±2.0 


+12.49 ±0.28 


HAT-P-23 














DS 


2008 May 19-2008 Sep 14 


5 


6000 ± 125 


4.5 ±0.25 


9.0 ±1.0 


-15.10 ±0.30 



^ The mean heliocentric RV of the target (in the lAU system). The error gives the rms of the individual velocity measures 
for the target with the given instrument. 



TABLE 3 

Summary of high-resolution/high-SN 

SPECTROSCOPIC observations 



Instrument 



Date(s) 



Number of 
RV obs. 



HAT-P-20 






Keck/HIRES 


2009 Apr-2009 Dec 


10 


HAT-P-21 






Keck/HIRES 


2009 May-2010 Feb 


15 


HAT-P-22 






Keck/HIRES 


2009 Apr-2009 Dec 


12 


HAT-P-23 






Keck/HIRES 


2008 Jun-2009 Dec 


13 



scription given by 'Vaughan. Preston fc Wilsonl ()1978D . 
and as described in .Hartman ct al . (2001. Note that 
o ur relative S index has not been cal ibrated to the scale 
of lVaughan. Preston k. Wilsonl (|1978[ ). We do not detect 
any significant variation of the index correlated with or- 
bital phase; such a correlation might have indicated that 
the RV variations could be due to stellar activity, casting 
doubt on the planetary nature of the candidate. There 
is no sign of emission in the cores of the Ca II H and K 
lines in any of our spectra, from which we conclude that 
all of the targets have low chromospheric activity levels. 

2.4. Photometric follow-up observations 

In order to permit a more accurate modeling of 
the light curves, we conducted additional photomet- 
ric observations with the KeplerCam CCD camera on 
the FLWO 1.2 m telescope for each star, and with the 
Faulkes North Telescope (FTN) of the Las Cumbres Ob- 
servatory Global Network (LCOGT) at Hawaii for HAT- 
P-21 only. The observations for each target are summa- 
rized in Table [H 

The reduction of these images, including basic cali- 
bration, astrometry, and apertu re photome try, was per- 
formed as described by Bakos et al.l (|201Cll ). We found 
that the aperture photometry for HAT-P-20 was signif- 
icantly affected by the close-by neighbor star 2MASS 
07273995+2420118 with Ai = 1.1 magnitude difference 
at 6.86" separation (Figure [5]). Thus, we performed im- 
age subtraction on the FLWO 1.2 m images with the 
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Fig. 7.— Keck/HIRES observations of HAT-P-21. The panels 
are as in Figu re [61 The parameters used in the best-fit model are 
given in Table [141 the RV jitter was 26.4ms~^, and the fit rms was 
26.6ms~'^. Observations shown twice are represented with open 
symbols. 

same toolset used for the HATNet reductions, but ap- 
plied a discrete kernel with half-size of 5 pixels and no 
spatial variations. Indeed, for this stellar configuration, 
the image subtraction results proved to be superior to 
the aperture photometry. For all of the follow-up light 
curves, we performed EPD and TFA to remove trends 
simultaneously with the lig ht curve modeling (for more 
details, see Section H and iBakos etall [2010) . The fi- 
nal time series, together with our best-fit transit light 
curve model, are shown in the top portion of Figures 1101 
through [131 ttie individual measurements are reported in 
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TABLE 4 

Relative radial velocities, bisector spans, and activity index 
measurements of hat-p-20. 



BJDutc 


RV'' 


o-Rv'' 


BS 


CBS 




CS 


(2,454,000+) 


(vn s 


(vn s 










934.84149 . 


-1139.51 


2.94 










954.82485 .. 






-8.00 


3.22 


2.051 


0.075 


954.83175 .. 


-942.88 


2.24 


-18.65 


8.02 


2.900 


0.035 


956.80481 .. 


1092.32 


2.28 


-42.21 


10.81 


3.305 


0.018 


983.73851 .. 


-1151.02 


2.39 


18.53 


4.51 


2.853 


0.016 


984.74064 .. 


328.09 


2.14 


6.85 


6.11 


2.879 


0.016 


986.74296 .. 


-1225.77 


2.34 


15.28 


7.78 


2.588 


0.029 


1107.15226 . 


-870.02 


4.89 


10.00 


17.65 


1.352 


0.043 


1109.05752 . 


1233.55 


3.76 


-1.23 


7.66 


2.783 


0.018 


1191.10274 . 


-1147.72 


2.07 


14.55 


6.16 


3.412 


0.016 


1192.11754 . 


1100.72 


2.12 


4.88 


2.83 


3.346 


0.015 


Note. — Note that for the iodine-free template exposures we do not 


measure 



the RV but do measure the BS and S index. Such template exposures can be 
distinguished by the missing RV value. 

^ Barycentric Julian dates throughout the paper are calculated from Coordi- 
nated Universal Time (UTC). 

^ The zero-point of these velocities is arbitrary. An overall offset 7rci fitted to 
these velocities in Section [33] has not been subtracted. 

^ Internal errors excluding the component of astrophysical jitter considered in 
Section |3.3l 

^ Relative chromospheric activity index, not calibrated to the scale of 
IVaughan. Preston & WilsonI i ll 9781) . 



TABLE 5 

Relative radial velocities, bisector spans, and activity index 
measurements of hat-p-21. 



BJDutc 


RV=' 


CRV 


b 


BS 


0"BS 




s= 




(2,454,000+) 


(m s 


^) 


(m s~ 




(m s 




(m s 










954.91548 . 










3 


67 


5 


10 





819 


0.002 


954.92307 .. 


221 


51 


2 


88 


-9 


60 


13 


62 





829 


0.003 


955.86532 .. 


-274 


01 


2 


48 


-8 


08 


5 


81 





823 


0.003 


955.98849 .. 


-314 


29 


2 


76 


7 


51 


12 


62 





855 


0.005 


957.00006 .. 


-378 


18 


3 


00 


2 


88 


8 


36 





866 


0.009 


963.95573 .. 


-197 


70 


2 


68 


-12 


46 


5 


48 





789 


0.003 


983.85638 .. 


212 


70 


2 


72 


-4 


61 


3 


33 





815 


0.003 


985.87988 .. 


-406 


23 


2 


46 


5 


88 


4 


95 





809 


0.003 


986.85645 .. 


509 


76 


3 


01 


-0 


70 


7 


24 





803 


0.005 


987.86480 .. 


270 


66 


2 


92 


-6 


96 


4 


20 





810 


0.004 


1174.14589 . 


-145 


85 


3 


13 





01 


6 


54 





833 


0.005 


1191.14119 . 


-391 


67 


2 


75 


5 


76 


5 


16 





832 


0.005 


1192.12302 . 


-457 


57 


2 


67 


7 


04 


3 


15 





809 


0.005 


1193.11001 . 


496 


07 


2 


58 


2 


50 


3 


85 





794 


0.005 


1197.15057 . 


459 


61 


2 


59 


9 


69 


4 


91 





777 


0.005 


1229.12322 . 


-469 


62 


2 


81 


-2 


53 


5 


81 





821 


0.005 



Note. — Notes for this table arc identical to that of Table [4] 



Tables [HHni for HAT-P-20 through HAT-P-23, respec- 
tively. 

3. ANALYSIS 

3.1. Properties of the parent stars 

Fundamental parameters for each of the host stars, 
including the mass (M^,) and radius (Ri,), which are 
needed to infer the planetary properties, depend strongly 
on other stellar quantities that can be derived spec- 
troscopically. For this we have relied on our template 
spectra obtained with the Keck/HIRES instrument, and 
the an alysis package known as S pectroscopy Made Easy 
(SME: iValenti fc Piskunovlll996D. along with the atomic 
line database o f IValenti fc Fischeil ([20051) . For each star, 
SME yielded the following initial values and uncertain- 
ties (which we have conservatively increased to include 



our estimates of the systematic errors): 

• HAT-P-20 - effective temperature Tcff* = 4626 ± 
104 K, stellar surface gravity logg* = 4.80 ± 
0.06 (cgs), metallicity [Fe/H] = -h0.31 ± 0.08 dex, 
and projected rotational velocity wsini = 3.1 ± 
0.5kms^^. 

• HAT-P-21 - effective temperature Toff* = 5701 ± 
144 K, stellar surface gravity logg^ = 4.48 ± 
0.14 (cgs), metallicity [Fe/H] = -1-0.04 ± 0.1 dex, 
and projected rotational velocity wsin?' = 4.1 ± 
0.5kms~^. 

• HAT-P-22 - effective temperature Toff* = 5338 ± 
88 K, stellar surface gravity logg^, = 4.52 ± 
0.14 (cgs), metallicity [Fe/H] = -^0.26 ± 0.08 dex, 
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TABLE 6 

Relative radial velocities, bisector spans, and activity index 
measurements of hat-p-22. 



BJDutc 


RV^ 


CRV 


b 


BS 


0"BS 








(2,454,000+) 


(ms^-*-) 


(m 




(m s 




(m s 










928.94934 . 


187.55 


1 


12 


13 


71 


4 


73 





557 


0.002 


954.88137 . . 


284.38 


1 


12 


7 


07 


3 


69 





556 


0.002 


954.88605 .. 








-2 


17 


3 


34 





554 


0.001 


955.79962 .. 


77.44 


1 


46 


6 


00 


13 


88 





533 


0.008 


956.96383 .. 


-273.80 


1 


30 


4 


41 


10 


92 





531 


0.005 


963.92157 .. 


-14.63 


1 


33 


-0 


17 


6 


57 





570 


0.002 


983.82710 .. 


282.58 





98 


2 


20 


3 


79 





572 


0.003 


985.81064 .. 


-298.55 


1 


14 


-2 


81 


2 


96 





560 


0.002 


986.84124 .. 


239.74 


1 


39 


-7 


48 


6 


58 





553 


0.003 


988.83025 .. 


-300.15 


1 


31 


-16 


80 


8 


24 





555 


0.002 


1191.14866 . 


-326.24 


1 


32 


-0 


33 


5 


00 





843 


0.003 


1193.10170 . 


255.62 


1 


30 


-8 


71 


6 


79 





835 


0.003 


1193.93435 . 


-192.43 


1 


38 


5 


08 


5 


27 





843 


0.004 



Note. — Notes for this tabic arc identical to that of Tabic [4] 



TABLE 7 

Relative radial velocities, bisector spans, and activity index 
measurements of hat-p-23. 



BJDutc 


RV^ 


crRV 


b 


BS 


0"BS 










(2,454,000+) 


(ms-l) 


(ms~ 




(m s 


') 


(m s 












638.09243 . 








-4 


89 


4 


59 


1 


170 





020 


638.10601 .. 


-137.90 


4 


49 


-0 


45 


5 


22 





797 





006 


674.91404 .. 


374.64 


5 


15 


-10 


71 


3 


99 





794 





006 


723.78719 .. 


-209.16 


5 


55 


-1 


68 


13 


38 





738 





006 


725.85309 .. 


329.57 


6 


18 


12 


86 


6 


55 





749 





006 


726.87976 .. 


277.98 


4 


20 


-26 


82 


6 


17 





770 





005 


727.76247 .. 


-210.01 


4 


81 


-35 


11 


5 


18 





811 





006 


727.89619 .. 


-19.91 


4 


77 


-30 


29 


5 


05 





795 





005 


777.82189 .. 


214.73 


4 


54 


-12 


64 


4 


43 





754 





005 


810.72960 .. 


290.22 


5 


63 


-13 


57 


5 


59 





726 





006 


955.04676 .. 


388.82 


4 


74 


-15 


76 


3 


57 





761 





006 


985.99494 .. 


-369.93 


4 


42 


-19 


69 


5 


06 





718 





006 


1192.70812 . 


348.01 


4 


40 


-14 


40 


4 


53 





726 





005 


1193.71080 . 


83.68 


4 


57 


-20 


67 


4 


23 





825 





008 



Note. — Notes for this table arc identical to that of Table [4] 



TABLE 8 

High-precision differential photometry of HAT-P-20. 



TABLE 9 

High-precision differential photometry of HAT-P-21. 



BJDutc 


Mag^ 


0"Mag 


Mag(orig)'^ Filter 


BJDutc 


Mag'' 


0"Mag 


Mag(orig)'^ Filter 


(2,400,000+) 








(2,400,000+) 






54902.59749 


-0.00169 


0.00075 


10.66286 i 


54942.76369 


0.01133 


0.00088 


10.43700 i 


54902.59799 


-0.00139 


0.00075 


10.66200 i 


54942.76432 


0.01370 


0.00088 


10.43950 i 


54902.59848 


-0.00469 


0.00075 


10.65893 i 


54942.76495 


0.00989 


0.00088 


10.43550 i 


54902.59919 


0.00281 


0.00075 


10.66596 i 


54942.76574 


0.00891 


0.00088 


10.43490 i 


54902.59970 


-0.00108 


0.00075 


10.66316 i 


54942.76637 


0.00974 


0.00088 


10.43570 i 


54902.60020 


0.00141 


0.00075 


10.66524 i 


54942.76698 


0.01073 


0.00088 


10.43650 i 


54902.60089 


0.00440 


0.00075 


10.66726 i 


54942.76779 


0.00961 


0.00088 


10.43480 i 


54902.60139 


0.00040 


0.00075 


10.66424 i 


54942.76839 


0.00938 


0.00088 


10.43450 i 


54902.60190 


-0.00064 


0.00075 


10.66234 i 


54942.76902 


0.01249 


0.00088 


10.43800 i 


54902.60259 


-0.00009 


0.00075 


10.66388 i 


54942.76980 


0.01011 


0.00088 


10.43520 i 



Note. — This table is available in a machine-readable form in the 
online journal. A portion is shown here for guidance regarding its form 
and content. 

^ The out-of-transit level has been subtracted. These magnitudes have 
been subjected to the EPD and TFA procedures, carried out simulta- 
neously with the transit fit. 

^ Raw magnitude values without application of the EPD and TFA 
procedures. 



Notes for this table are identical to those of Table [sl 



Note. 



HAT-P-23 - effective temperature Tcff* = 5987 ± 
120 K, stellar surface gravity logg* — 4.48 ± 
0.12 (cgs), metallicity [Fe/H] = -^0.18 ± 0.1 dex, 
and projected rotational velocity wsini ~ 8.0 ± 
0.5kins~^. 



and projected rotational velocity wsini = 1.5 ± 
0.5 kms^^ 



In principle the effective temperatiue and metallicity, 
along with the surface gravity taken as a luminosity in- 
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Fig. 8.— Keck/HIRES observations of HAT-P-22. The panels 
are as in Figu re \§ \ The parameters used in the best-fit model are 
given in Table [T4l the RV jitter was 9.9ms~^, and the fit rms was 
10.0 ms^^. Observations shown twice are represented with open 
symbols. 
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40 L 
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Phase with respect to 

Fig. 9.— Keck/HIRES observations of HAT-P-23. The panels 
are as in Figu re [61 The parameters used in the best-fit model are 
given in Table [Ml and the RV jitter was 34.7 ms~^. Observations 
shown twice are represented with open symbols. 



TABLE 10 

High-precision differential photometr.y of IfAT-P-22. 



2009 Uar 1 1 

'l/o'.Jrt.^.'rVy'; 



BJDuTC 
(2,400,000-1-) 


Mag^ 


0"Mag 


Mag(orig)'^ Filter 


54891.60570 


0.00028 


0.00069 


8.78803 


z 


54891.60603 


-0.00100 


0.00069 


8.78658 


z 


54891.60636 


-0.00189 


0.00069 


8.78541 


z 


54891.60667 


-0.00138 


0.00069 


8.78514 


z 


54891.60702 


-0.00237 


0.00069 


8.78534 


z 


54891.60735 


0.00030 


0.00069 


8.78923 


z 


54891.60767 


0.00404 


0.00069 


8.79234 


z 


54891.60801 


-0.00170 


0.00069 


8.78696 


z 


54891.60831 


0.00257 


0.00069 


8.79009 


z 


54891.60865 


-0.00145 


0.00069 


8.78653 


z 


Note. — Notes for this table are 


identical to those of Table [sl 





dicator, could be used as constraints to infer the stellar 
mass and radius by comparison with stellar evolution 
models. However, the effect of logg* on the spectral line 
shapes is rather subtle, and as a result it is typically 
difficult to determine accurately, so that it is a rather 
poor luminosity indicator in practice. Unfortunately a 
trigonometric parallax is not available for any of the host 
stars, since they were not included among th e targets of 
the Hipparcos mission (|Perryman et al1ll997l ). For plan- 
etary transits, another possible constraint is provided by 
the a/Ri, normalized semi-major axis, which is closely 
related to the mean stellar density. The quantity 
a/Ri, can be derived directly from the combination of 





-0.1 -0.05 0.05 0.1 

Time from transit center (days) 

Fig. 10. — Unbinned transit light curves for HAT-P-20, acquired 
with KeplerCam at the FLWO 1.2 m telescope. The li ght c urves 
have been EPD and TEA processed, as described in S 13.31 The 
dates of the events are indicated. Curves after the first are dis- 
placed vertically for clar ity. Our best fit from the global modeling 
described in Section l3.3l is shown by the solid lines. Residuals from 
the fits are displayed at the bottom, in the same order as the top 
curves. The error bars represent the photon and background shot 
noise, plus the readout noise. 

the transit light curves (jSozzetti et al.|[2007l) and the RV 
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2009. Apr 20 
2010 Feb 15 



• 2010 Fab 19 
• »• • • * • 




-0.05 0.05 

Time from transit center (days) 

Fig. 11. — Similar to Figure llOl here we show the follow-up light 
curves for HAT-P-21. 



£ 0.02 



2009 Feb 28 





-0.1 -0.05 0.05 0.1 

Time from transit center (days) 

Fig. 12. — Similar to Figure fTOl here we show the follow-up light 
curves for HAT-P-22. 



data (required for eccentric cases, see Section ESI) • This, 
in turn, allows us to improve on the determination of the 
spectroscopic parameters by supplying an indirect con- 
straint on the weakly determined spectroscopic value of 
log (7* that removes degeneracies. We take this approach 
here, as described below. The validity of our assumption, 
namely that the adequate physical model describing our 
data is a planetary transit (as opposed to a blend), is 
shown later in Section |3^ 




2008 Jun 14 



2008 Sep 08 



2008 Sep 13 



2009 Apr 19 



3D09 Jui 13 



irjf^lf^-'ii^isi^^f^i^Jfg^'ifs^isf^i- 

I 

-^xg^g^^igiigfli^iei^^iiif^^if 

-it'— •iitH^ i-rjS5nS«il^ft!^j<^:sJl^ 
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-0.1 -0.05 0.05 0.1 

Time from transit center (days) 

Fig. 13. — Similar to Figure fTOl here we show the follow-up light 
curves for HAT-P-23. 



TABLE 11 

High-precision differential photometry of HAT-P-23. 



BJDuTC 


Mag^ 


0"Mag 


Mag(orig)'' Filter 


(2,400,000-1-) 






54632.73784 


0.02143 


0.00111 


10.98360 i 


54632.73870 


0.01830 


0.00116 


10.98170 i 


54632.73954 


0.01503 


0.00115 


10.97810 i 


54632.74039 


0.01859 


0.00115 


10.98340 i 


54632.74125 


0.01123 


0.00111 


10.97530 i 


54632.74208 


0.02006 


0.00109 


10.98390 i 


54632.74295 


0.01831 


0.00111 


10.98220 i 


54632.74381 


0.01766 


0.00114 


10.97840 i 


54632.74637 


0.02025 


0.00117 


10.99540 i 


54632.74722 


0.01556 


0.00112 


10.97990 i 


Note. — Notes for this table 


are identical to those of Table Is] 



For each system, our initial values of Toff*, log 5*, 
and [Fe/H] were used to determine auxiliary quantities 
needed in the global modeling of the follow-up photome- 
try and radial velocities (specifically, the limb-darkening 
coefficients). This modeling, the details of which are 
described in Section [331 uses a Monte Carlo approach 
to deliver the numerical probability distribution of a/Ri, 
and other fit ted v ariable s. For further details we refer 
the reader to iPall ()2009bD . When combining a/i?* (used 
as a proxy for luminosity) with assumed Gaussian dis- 
tributions for Teffi, and [Fe/H] based on the SME deter- 
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minations, a comparison with stellar evolution models 
allows the probability distributions of other stellar prop- 
erties to be inferred, including logg*. Here we use the 
stellar evolution calcu lations from the Yonsei-Yale group 
(YY; lYi et al.l l2001| ) for all planets presented in this 
work. The comparison against the model isochrones was 
carried out for each of 10,000 Monte Carlo trial sets for 
HAT-P-21, HAT-P-22, and HAT-P-23, and for 20,000 
Monte Carlo trial sets for HAT-P-20 (see Section [53)1 . 
Parameter combinations corresponding to unphysical lo- 
cations in the H-R diagram (26% of the trials for HAT- 
P-20, and less than 1% of the trials for the other objects) 
were ignored, and replaced with another randomly drawn 
parameter set. For each system we carried out a second 
SME iteration in which we adopted the value of \ogg^, 
so determined and held it fixed in a new SME analysis 
(coupled with a new global modeling of the RV and light 
curves), adjusting only Tch*, [Fe/H], and i;sini. This 
gave: 

• HAT-P-20: \ogg^ = 4.64 ± 0.06, T^s* = 4595 ± 
80 K, [Fe/H] = -hO.35 ± 0.08, and vsini = 2.1 ± 
O.Skms-i. 

• HAT-P-21: log 5* = 4.31 ± 0.06, Toff* = 5588 ± 
80 K, [Fe/H] = +0.01 ± 0.08, and vsini = 3.5 ± 
0.5 kms^^. 

• HAT-P-22: \ogg^ 4.37 ± 0.06, Teff* = 5302 ± 
80 K, [Fc/H] = +0.24 ± 0.08, and vsini = 0.5 ± 
0.5kms-i. 

• HAT-P-23: \ogg^ = 4.33 ± 0.06, Toff* = 5905 ± 
80 K, [Fe/H] = +0.15 ± 0.04, and vsini = 8.1 ± 
0.5kms-^ 

In each case the conservative uncertainties for Toff* and 
[Fe/H] have been increased by a factor of two over their 
formal values, as before. For each system, a further itera- 
tion did not change log significantly, so we adopted the 
values stated above as the final atmospheric properties 
of the stars. They are collected in Table [T2] 

With the adopted spectroscopic parameters the model 
isochrones yield the stellar mass and radius, and other 
properties. These are listed for each of the systems in 
Table [HI According to these models HAT-P-20 is a 
dwarf star with an estimated age of 6.7^3 8 Gyr, HAT- 
P-21 is a slightly evolved star with an estimated age 
of 10.2 ± 2.5 Gyr, HAT-P-22 is a slightly evolved star 
with an estimated age of 12.4 ± 2.6 Gyr, and HAT-P- 
23 is a slightly evolved star with an estimated age of 
4.0 ±1.0 Gyr. The inferred location of each star in a dia- 
gram of a/ Ri, versus Toff*, analogous to the classical H-R 
diagram, is shown in Figure [TH In all cases the stellar 
properties and their la and 2a confidence ellipsoides are 
displayed against the backdrop of model isochrones for 
a range of ages, and the appropriate stellar metallicity. 
For comparison, the locations implied by the initial SME 
results are also shown (in each case with a triangle). 

The stellar evolution modeling provides color indices 
(see Table fT2|) that may be compared against the mea- 
sured values as a sanity check. For each star, the 
best available measurements are t he near-infrared magn i- 
tudes from the 2MASS Catalogue (jSkrutskie et al.ll2006[ ). 
which are given in Table [12] These are converted to the 



photometric system of the models (ESQ system) using 
the transformations bv iCarpenteii ()2001[ ). The resulting 
color indices are also shown in Table [T2] for HAT-P-20 
through HAT-P-23, respectively. Indeed, the colors from 
the stellar evolution models and from the observations 
agree for all of the host stars within 2-a. The distance 
to each object may be computed from the absolute K 
magnitude from the models and the 2MASS Kg mag- 
nitudes, which has the advantage of being less affected 
by extinction than optical magnitudes. The results are 
given in Table [121 where in all cases the uncertainty ex- 
cludes possible systematics in the model isochrones that 
are difficult to quantify. 

3.2. Excluding blend scenarios 

Our initial spectroscopic analyses discussed in 
Section [2?2] and Section [2731 rule out the most obvious as- 
trophysical false positive scenarios. However, more sub- 
tle phenomena such as blends (contamination by an un- 
resolved eclipsing binary, whether in the background or 
associated with the target) can still mimic both the pho- 
tometric and spectroscopic signatures we see. In the fol- 
lowing section we investigate whether such scenarios may 
have caused the observed photometric and spectroscopic 
features. 

3.2.1. Spectral line-bisector analysis 

Following [Torres et al.l (|2007l ). we explored the possi- 
bility that the measured radial velocities are not real, 
but are instead caused by distortions in the spectral 
line profiles due to contamination from a nearby unre- 
solved eclipsing binary. A bisector span (BS) analysis 
for each system based on the K eck spectra was done as 
described in §5 of Bak os et al.l ([2007'). In general, none 
of the Kcck/HIRES spectra suffer significant sky con- 
tamination. Nevertheless, we calculate d the Sky Con- 
tamin ation Factors (SCF) as described in lHartman et al] 
(|2009f) . and corrected for the minor correlation between 
SCF and BS. The results are exhibited in Figure [TSl 
where we show the SCF-BS and RV-BSscf (BS after 
SCF correction) plots for each planetary system. We 
also calculated the Spearman rank-order correlation co- 
efficients (denoted as Rs for the RV vs. BS quantities 
and the false alarm probabilities (see Table [T3]). There is 
no correlation for HAT-P-21, HAT-P-22 and HAT-P-23, 
and thus the interpretation of these systems as transit- 
ing planets is clear. There is an anti-correlation present 
for HAT-P-20, which is strengthened when the SCF cor- 
rection is applied. A plausible explanation for this is 
that the neighboring star at 6" separation (see Figure [5]) 
is bleeding into the slit, even though we were careful 
during the observations to keep the slit centered on the 
main target, and adjusted the slit orientation to be per- 
pendicular to the direction to the neighbor. We simu- 
lated this scenario, and calculated the expected BS as a 
function of RV due to the neighbor, assuming that the 
two stars have the same systemic velocity and the see- 
ing is 1". Indeed, we get a slight anti-correlation from 
this simulation, and the range of magnitude in the BS 
variation is consistent with the observations. It is also 
possible that some of the anti-correlation is due to the 
fact that the slit was not in vertical angle for many of 
the observations. The non-vertical slit mode may result 
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Fig. 14. — Upper left: Model isochrones from lYi et al.l l|2001l ) for the measured metallicity of HAT-P-20, [Fe/H]= +0.35, and ages of 
3.0, 4.0, 5.0, 5.5, 6.0, 6.5, 7.0, 7.5, 8.0, 8.5, 9.0, 10.0 and 11.0 Gyr (left to right). The adopted values of T^s^, and a/ij* are shown together 
with their Icr and 2ct confidence ellipsoids. The initial values of T^ff-t. and a/Ri, from the first SME and light curve analyses are represented 
with a triangle. Upper right: Same as upper left, here we show the results for HAT-P-21, with [Fe/H]= +0.01, and ages of 5.0, 6.0, 7.0, 
8.0, 9.0, 10.0, 11.0, 12.0, and 13.0 Gyr (left to right). Lower left: Same as upper left, here we show the results for HAT-P-22, with [Fe/H] = 
+0.24, and ages of 5.0, 6.0, 7.0, 8.0, 9.0, 10., 11.0, 12., 13.0, and 14.0 Gyr (left to right). Lower right: Same as upper left, here we show the 
results for HAT-P-23, with [Fe/H]= +0.15, and ages of 0.2, 0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5 and 5.0 Gyr (left to right). 



in wavelength-dependent slit losses due to atmospheric 
dispersion, and this could bring in a correlation with the 
sky background, and change the shape of the spectral 
lines. 



3.3. Global modeling of the data 

This section describes the procedure we followed for 
each system to model the HATNet photometry, the 
follow-up photometry, and the radial velocities simulta- 
neously. Our model for the follow-up light curves used 
analytic formulae based on iMandel fc A gol' (2002") for 
the eclipse of a star by a planet, with limb darkening 
being prescribed by a quadratic law. The limb dark- 
ening coefficients for the Sloan g band, Sloan i band, 
and Sl o an z band were interpolated from the tables by 
iClaretl (|2004D for the spectroscopic parameters of each 
star as determined from the SME analysis (Section 13. ip . 
The transit shape was parametrized by the normalized 
planetary radius p = Rp/Ri,, the square of the impact 
parameter 6^, and the reciprocal of the half duration of 
the transit C/R*. We chose these parameters because of 
their simple geometric meaning s and the fact that these 
show negligible correlations (see iBakos et al.|[20lol ). The 
relation between C/R* and the quantity a/Rt, used in 
Section [3TT| is given by 



a/R^, = P/27r(C/i?*)\/l - 6Vl - e^/{l + esmuj) (1) 



(see, e.g-. lTinglev fc SackettI 120051 ). Note the subtle de- 
pendency of a/Ri, on the k = e cos w and h = e sin uj 
Lagrangian orbital parameters that are typically derived 
from the RV data {lu is the longitude of periastron). This 
dependency is often ignored in the literature, and a/i?* 
is quoted as a "pure" light curve parameter. Of course, 
if high quality secondary eclipse observations are avail- 
able that determine both the location and duration of 
the occultation, then k and h can be determined with- 
out RV data. Our model for the HATNet data was the 
simplified "P1P3" version of the IMandel fc Agol (pOOl 
analytic functions (an expansion in terms of Legendre 
polyn omials), for the reasons described in Bakos et al] 
(2010'). Following the formahsm presented bv[Pall ([2009), 
the RVs were fitted with an eccentric Keplerian model 
parametrized by the semi- amplitude K and Lagrangian 
elements k and h. Note that we allowed for an eccentric 
orbit for all planets, even if the results were consistent 
with a circular orbit. There are several reasons for this: 
i) many of the close-in hot Jupiters show eccentric or- 
bits, thus the assumption of fixing e = has no physical 
justification (while this has been customary in early dis- 
coveries relying on very few data-points) ii) the error-bars 
on various other derived quantities (including a/i?*) are 
more realistic with the inclusion of eccentricity, and iii) 
non-zero eccentricities can be very important in proper 
interpretation of these systems. 
We assumed that there is a strict periodicity in the 
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TABLE 12 

Stellar parameters for HAT-P-20-HAT-P-23 



Parameter HAT-P-20 HAT-P-21 HAT-P-22 HAT-P-23 Source 



Spectroscopic properties 

Tcff* (K) 4595 ± 80 

[Fe/H] +0.35 ± 0.08 

V sini (km s~^)... 2.1±0.5 

VmiLc (kms~^) 2.21 

Vmic (kms-i) .... 0.85 

7RV (kms-i) -18.81 ±0.68 

Photometric properties 

y (mag) 11.339 

V-Ic (mag) 1.50 ±0.12 

J (mag) 9.276 ±0.022 

H (mag) 8.743 ±0.021 

Ks (mag) 8.601 ± 0.019 

J - K (mag,ESO) 0.715 ± 0.033 

Derived properties 

(Mq) 0.756 ±0.028 

i?* (Rq) 0.694 ±0.021 

logg* (cgs) 4.63 ±0.02 

L^, (Lq) 0.19 ±0.02 

Mv (mag) 7.07 ±0.17 

Mk (mag,ESO) . . 4.42 ± 0.09 
J- K (mag.ESO) 0.66 ± 0.02 

Age (Gyr) 6.7^11 

Distance (pc) 70 ± 3 



5588 ± 80 
+0.01 ±0.08 
3.5 ± 0.5 
3.74 
0.85 
-53.19 ±0.09 



11.685 
0.654 ±0.097 
10.503 ±0.022 
10.154 ±0.019 
10.111 ±0.018 
0.419 ±0.031 



0.947 ± 0.042 
1.105 ± 0.083 
4.33 ± 0.06 
1 Q6+°-^° 

^•'^"-0.16 

4.80 ±0.19 
3.12 ± 0.16 
0.44 ± 0.02 
10.2 ± 2.5 
254 ± 19 



5302 ± 80 
+0.24 + 0.08 
0.5 ± 0.5 
3.30 
0.85 
+12.49 + 0.28 



9.732 
0.992 + 0.066 
8.293 + 0.023 
7.935 + 0.029 
7.837 + 0.021 
0.486 + 0.034 



0.916 + 0.035 
1.040 + 0.044 
4.36 + 0.04 
0.77 + 0.09 
5.22 + 0.14 
3.30 + 0.10 
0.50 + 0.02 
12.4 + 2.6 
82 + 3 



5905 + 80 SME'' 

fO.15 + 0.04 SME 

8.1 + 0.5 SME 

4.22 SME 

0.85 SME 

-15.10 + 0.30 DS 



12.432 
0.68 + 0.12 
11.103 + 0.022 
10.846 + 0.022 
10.791 + 0.020 
0.335 + 0.032 



1.130 + 0.035 
1.203 + 0.074 
4.33 + 0.05 
1.58 + 0.23 
4.31 + 0.16 
2.86 + 0.14 
0.37 + 0.01 
4.0 + 1.0 
393 + 25 



TASS 

TASS 

2MASS 

2MASS 

2MASS 

2MASS 



YY+a/R. 
YY+a/R. 
YY+a/R, 
YY+a/R, 
YY+a/R. 
YY+a/R. 
YY+a/R, 
YY+a/R, 
YY+a/R 



-+SME b 

-+SME 

-+SME 

-+SME 

-+SME 

-+SME 

-+SME 

-+SME 

.+SME 



^ SME — "Spectroscopy Made Easy" package for the analysis of high-resolution spectra jValenti & Piskunoviri996D . 
These parameters rely primarily on SME, but have a small dependence also on the iterative analysis incorporating the 
isochrone search and global modeling of the data, as described in the text. 

^ YY+a/H*+SME ^ Based on the YY isochrones \Y\ et al.||200lD . a/i?* as a luminosity indicator, and the SME results. 



TABLE 13 
Summary of RV vs. BS correlations. 



Name 




FAPi*^ 




FAP2'* 


HAT-P-20 


-0.73 


2.46% 


-0.92 


0.05% 


HAT-P-21 


-0.24 


39% 


-0.20 


47% 


HAT-P-22 


0.33 


30% 


0.27 


37% 


HAT-P-23 


0.27 


37% 


0.27 


37% 



^ The Spearman correlation coefficient between 
tlie bisector (BS) variations and the radial veloc- 
ities (RV). 



False alarm probability for Rsl- 

The Spearman correlation coefficient between 

BS corrected for the sky contamination factor 

(SCF) and the RVs. 

False alarm probability for i?,s2. 

individual transit times. For each system we assigned 
the transit number Ntr = to a complete follow-up 
light curve. For HAT-P-20b this was the light curve 
gathered on 2009 Oct 21, for HAT-P-21b: 2010 Feb 19, 
HAT-P-22b: 2009 Feb 28, and HAT-P-23b: 2008 Sep 13. 
The adjustable parameters in the fit that determine the 
ephemeris were chosen to be the time of the first transit 
center observed with HATNet (Tc -252, 7c -286, Tc -135, 
and Tc _3i2 for HAT-P-20b through HAT-'P-23b respec- 
tively) and that of the last transit center observed with 
the FLWO 1.2m telescope (Tco, Td, 19, and 250 
for HAT-P-20b through HAT-P-23b respectively). We 
used these as opposed to period and reference epoch in 
order to minim ize correlations between parameters (see 
iPal et aLll2Q08D . Times of mid-transit for intermediate 
events were interpolated using these two epochs and the 



corresponding transit number of each event, Ntr- The 
eight main parameters describing the physical model for 
each system were thus the first and last transit center 
times, Rp/Ri,, 6^, K, k = ecosw, and h = esinoj. 

For HAT-P-20b, HAT-P-22b, and HAT-P-23b three ad- 
ditional parameters were included (for each system) that 
have to do with the instrumental configuration (blend 
factor, out-of-transit magnitudes, gamma velocities; see 
later). For HAT-P-21b seven additional parameters were 
included, because it was observed in 3 different HATNet 
fields. These are the HATNet blend factor i?inst (one for 
each HATNet field), which accounts for possible dilution 
of the transit in the HATNet light curve from background 
stars due to the broad PSF (20" FWHM), the HATNet 
out-of-transit magnitude Mq, HATNet (one for each HAT- 
Net field), and the relative zero-point 7rei of the Keck 
RVs. 

We extended our physical model with an instrumental 
model that describes brightness variations caused by sys- 
tematic errors in the measurements. This was done in a 
similar fashion to the analysis presented bv iBakos et al.l 
(|2010D . The HATNet photometry has already been EPD- 
and TFA-corrected before the global modeling, so we 
only considered corrections for systematics in the follow- 
up light curves. We chose the "ELTG" method, i.e., 
EPD was performed in "local" mode with EPD coeffi- 
cients defined for each night, and TFA was performed in 
"global" mode using the same set of stars and TFA coef- 
ficients for all nights. The five EPD parameters were the 
hour angle (representing a monotonic trend that changes 
linearly over time) , the square of the hour angle (refiect- 
ing elevation), and the stellar profile parameters (equiv- 
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Fig. 15. — Panels on the left show the bisector spans (BS) as a function of Sky Contamination Fa<;tor (SCF). Panels on the right exhibit 
the SCF-eorrected BS as a function of the radial velocities. The individual planets are labeled. 



alent to FWHM, elongation, and position angle of the 
image). The functional forms of the above parameters 



contained six coefficients, including the auxiliary out-of- 
transit magnitude of the individual events. For each sys- 
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tern the EPD parameters were independent for all nights, 
implying 12, 18, 12, and 36 additional coefficients in 
the global fit for HAT-P-20b through HAT-P-23b respec- 
tively. For the global TEA analysis we chose 20, 3, 10, 
and 20 template stars for HAT-P-20b through HAT-P- 
23b that had good quality measurements for all nights 
and on all frames, implying an additional 20, 3, 10, and 
20 parameters in the fit for each system. In all cases 
the total number of fitted parameters (43, 36, 33 and 67 
for HAT-P-20b through HAT-P-23b) was much smaller 
than the number of data points (755, 1172, 892 and 953, 
counting only RV measurements and follow-up photom- 
etry measurements). 

The joint fit was performed as described in lBakos et al.l 
(|2010f) . We minimized in the space of parameters by 
using a hybrid algorithm, co mbining the down hill sim- 
plex method (AMOEBA; see iPress et al.l 119921 ) with a 
classical linear least squares algorithm. Uncertainties for 
the parameters were derived by applying the Markov 
Chain Monte-Carlo method (MCMC, see iFb^ l2006l) . 
This provided the full a posteriori probability distribu- 
tions of all adjusted variables. The a priori distribu- 
tions of the parameters for these chains were chosen to 
be Gaussian, with eigenvalues and eigenvectors derived 
from the Eisher covariance matrix for the best-fit solu- 
tion. The Eisher covariance matrix was calc ulat e d ana - 
lytically using the partial derivatives given bv lPgl(l200a . 

EoUowing this procedure we obtained the a posteriori 
distributions for all fitted variables, and other quantities 
of interest such as a/R^,. As described in Section [XTl 
a/Ri, was used together with stellar evolution models to 
infer a value for log gi, that is significantly more accurate 
than the spectroscopic value. The improved estimate was 
in turn applied to a second iteration of the SME anal- 
ysis, as explained previously, in order to obtain better 
estimates of Tcff* and [Ee/H]. The global modeling was 
then repeated with updated limb-darkening coefficients 
based on those new spectroscopic determinations. The 
resulting geometric parameters pertaining to the light 
curves and velocity curves for each system are listed in 
Table [H 

Included in each table is the RV "jitter", which is a 
noise term that we added in quadrature to the internal 
errors for the RVs in order to achieve x^/dof = 1 from 
the RV data for the global fit. The jitter is a combination 
of assumed astrophysical noise intrinsic to the star, plus 
instrumental noise rising from uncorrected intstrumental 
effects (such as a template spectrum taken under sub- 
optimal conditions). 

The planetary parameters and their uncertainties were 
derived by combining the a posteriori distributions for 
the stellar, light curve, and RV parameters. In this way 
wc find masses and radii for each planet. These and 
other planetary parameters are listed at the bottom of 
Table [T31 and further discussed in § 21 

4. DISCUSSION 

4.1. HAT-P-20b 

IIAT-P-20b is a very massive {Mp = 7.246 ± 
0.187 Mj = 2302.9 ± 59.5 M®) and very compact 
(Rr> = 0.8 67 ± 0.033 i?j) hot Jupiter orbiting a K3 
ISkifj I2OO90 star. HAT-P-20b is the sixth most mas- 
sive, and second most dense transiting planet with 



Pp =13.78 ± 1.50gcm-3(see Eigure HH). The only 
planet (or brown dwarf) denser than HAT-P-20b is 
CoRoT-3b (jPeleuil et al.l [20081 ) with pp w 27gcm-3. 
Modeling IIAT-P-20b may be a challenge, as the old- 
est (4Gyr, i . e. yie lding the most compact planets) 
lEortnev et all (|2007D models with Mp = 2154 to- 
tal mass and 100 core-mass predict a much bigger 
radius (1.04i?j). The observed radius of 0.87 Mj would 
require a very high metal content. We note that the 
host star is one of the most metal rich stars that have 
a transiting planet ([Ee/H] = -1-0.35 ± 0.08). Curiously, 
IIAT-P-20b orbits a fairly late type star (K3), as com- 
pared to most of the massive hot Jupiters that orbit ~E5 
dwarfs. It is also different from the rest of the population 
in that the orbit is consistent with circular at the 3(7 level. 
The irradiation IIAT-P-20b receives is on e of the small- 
est, clearly making it a pL class exoplanet ()Eortnev et al.l 
I2OOI : (F) = (2.00±0.19)-10*^ergs-icm-2, comparable 
to the mean flux per orbit for another "heavy" planet 
HD 17156b on a 21 d period orbit. The only other mas- 
sive planet that receives less average flux (integrated over 
an orbit) is HD 80606b. HAT-P-20 is an extreme outlier 
in the Alp-Mi, plane; it is a relatively small mass star 
harboring a very massive planet. Another outlier (al- 
beit to a much lesser extent) with similar planetary ra- 
dius and stellar mass is WASP-lOb (IJohnson et al.l[200l 
[Christian et al.ir2009 ^. but this planet has less than half 
of the mass of HAT-P-20b (3.09 Mj). We also calculated 
the maximum mass of a stable moon for both the pro- 
grade and retrograde cases, and derived 0.128 M0 and 
8.31 Mq, respectively, i.e. HAT-P-20b can harbor a fairly 
massive moon. An 8.31 M^ retrograde moon would cause 
~ 10 s variations in the transit times, which is marginally 
detectable from the ground. 

HAT-P-20 has a close-by faint and red companion star 
at ~ 6.86" separation. Based on the Palomar sky sur- 
vey archival plates, we confirm that they form a close 
common-proper motion pair, thus it is very likely that 
the two stars are physically associated. The binary has 
appeared in the Washington Double Star co mpilation 
WDS ) as POU2795, and was discovered by iPourteaul 
19331 ). Eurthermore, based on the summary of obser- 
vations in the WDS, there is already a hint of orbital 
motion of the companion to HAT-P-20 over the last cen- 
tury. The position angle of the companion changed from 
PA=323° to PA=320° over the course of 89 years (be- 
tween 1909 and 1998), and it seems to be retrograde on 
the sky (clockwise). Thus, HAT-P-20 is yet another ex- 
ampl e of a massive planet in a binary system (jUdrv et al.l 
\20m . The binary companion makes this system ideal for 
high precision ground or space-based studies, as it pro- 
vides a natural comparison source, even though it has a 
later spectral type. 

4.2. HAT-P-21h 

With a mass of Mp = 4.063 ± 0.161 Mj, HAT-P-21b 
is the 11th most massive transiting planet. HAT-P-21b 
has a radius of Rp — 1.024 ± 0.092 mean density 

Pp ~ 4.68lg gg gcm~'^, aud orbits on a moderately eccen- 
tric orbit with e = 0.228 ±0.016, w = 30 9 ±3°. The tran- 
sits oc cur near apastron. As noted by iBuchhave et al.l 
(|2010[ ) , 4 Mj mass planets are very rare in the sample of 
currently known transiting exoplanets, and the only sib- 
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TABLE 14 

Orbital and planetary parameters for HAT-P-20b-HAT-P-23b 



Parameter HAT-P-20b 

Light curve parameters 

P (days) 2.875317 ± 0.000004 

Tc (BJDuTc) 2455080.92661 ± 0.00021 

Ti4 (days) 0.0770 ± 0.0008 

Ti2 = T34 (days) 0.0137 ± 0.0009 

a/_R* 11.17 ±0.29 

C/K* 31.32 ±0.22 

Rp/Ri, 0.1284 ±0.0016 

0.398tHlI 

b = acosi/R-i, O-BSlloojl 

i (deg) 86.8 ± 0.2 

Limb-darkening coefficients ^ 

ai (linear term) 0.4719 

bi (quadratic term) 0.2174 

«9 

h 

RV parameters 

K (ms-l) 1246.0 ±8.1 

kBy" 0.012 ±0.004 

hjiY" -0.007 ±0.008 

e 0.015 ±0.005 

uj (deg) 317 ± 130 

RV jitter (ms~i) 15.7 

Secondary eclipse parameters (derived) 

Ts (BJDuTc) 2455082.385 ± 0.007 

Ts 14 0.0764 ± 0.0010 

Ts,i2 0.0134 ± 0.0009 

Planetary parameters 

Mp (Mj) 7.246 ±0.187 

Rp (Rj) 0.867 ±0.033 

C{Mp,Rp) 0.50 

Pp (gcm-3) 13.78 ± 1.50 

loggp (cgs) 4.38 ±0.03 

a (AU) 0.0361 ± 0.0005 

Tcq (K) 970 ± 23 

e'^ 0.794 ±0.031 

Fper (10*ergs-i cm-^) f . . . , 2.O6 ± 0.20 

Fap (10*ergs-i cm-2) i 1.94 ±0.19 

(F) (10*ergs-lcm-2) f 2.00 ±0.19 



HAT-P-21b 



HAT-P-22b 



HAT-P-23b 



4.124481 ±0.000007 
2454996.41312 ± 0.00069 
0.1530 ±0.0027 
0.0184 ± 0.0029 
9.60 ±0.71 
14.81 ± 0.15 
0.0950 ± 0.0022 

i+0.089 
-0.118 



0.298"' 



546+° °'''' 
87.2 ± 0.7 



0.2976 
0.3131 



548.3 ± 14.2 
0.147 ±0.011 
-0.175 ± 0.019 
0.228 ±0.016 
309 ±3 
26.4 



2454998.865 ± 0.029 
0.1163 ±0.0068 
0.0117 ±0.0015 



4.063 ±0.161 
1.024 ± 0.092 
0.28 

4-b8_(, gg 

3.98 ±0.08 
0.0494 ± 0.0007 
1283 ± 50 
0.413 ±0.038 
10.0 ± 1.5 
3.96 ±0.66 
6.12 ±0.97 



3.212220 ± 0.000009 
2454930.22001 ± 0.00025 
0.1196 ± 0.0014 
0.0144 ±0.0013 
8.55 ± 0.35 
18.97 ±0.09 
0.1065 ± 0.0017 

r, 9-174-0.052 
"■^^^'-0.065 
U.4DD_Q Qg3 
-,+0.6 
-0.5 



86.9^ 



0.3587 
0.2831 
0.7080 
0.1165 



313.3 ±4.2 
-0.008 ± 0.008 
0.004 ± 0.016 
0.016 ± 0.009 
156 ± 66 
9.9 



2454931.809 ± 0.016 
0.1204 ± 0.0032 
0.0145 ± 0.0014 



2.147 ±0.061 
1.080 ±0.058 
0.36 

9 11+0.40 
^•-^-^--0.29 

3.66 ± 0.05 
0.0414 ± 0.0005 
1283 ± 32 
0.179 ± 0.010 
6.33 ± 0.67 
5.91 ± 0.60 
6.12 ± 0.62 



1.212884 ± 0.000002 
2454852.26464 ± 0.00018 
0.0908 ± 0.0007 
0.0105 ± 0.0007 
4.14 ±0.23 
24.90 ±0.12 
0.1169 ±0.0012 

n iQc:+0.053 

U.iUO_Q Q4g 

+0.070 



0-324_oigi 
85.1 ± 1.5 



0.2524 
0.3426 
0.5311 
0.2539 



368.5 ± 17.6 
-0.048 ± 0.023 
0.090 ±0.052 
0.106 ±0.044 
118 ±25 
34.7 



2454852.834 ± 0.018 
0.1064 ±0.0095 
0.0129 ±0.0020 



2.090 ±0.111 
1.368 ±0.090 
0.56 
1.01 ±0.18 
3.44 ±0.05 
0.0232 ± 0.0002 
2056 ± 66 
0.062 ± 0.004 
50.0 ± 11.4 
32.7 ±2.7 
40.3 ±5.3 



Tc: Reference epoch 
ration, time between f 

Values for a quadratii 
Lagrangian orbital pa 



af mid transit that minimizes tire correlati 
:st to last contact; T12 = T34: ingresa/egrc 
law, adopted from the tabulations by HJlan 
ameters derived from the global modeling, 



( tin 



! between f 



elati< 



effic 



ind prima 



zcording to the spectroscopic (SME) p 
ily determined by the RV data. 



JV^^ = -16. B.1D i£ 
i and fourth contact. 

ters listed in Table fT2l 



Iculated from UTC. T14: total tr; 



The Safronov nun 
Incoming flux per 



nt between the planetary mass Mp and r 
r is given by = ^ (V^so / V^rbf = (0/ Hp ) ( Mp /M, ) (s. 
it surface area. (F) is averaged over the orbit. 



lings of HAT-P-21b are HD 80606b (4.08 Afj; Naef et al.' 
[20011 ) and HAT-P-16b (4.19 M.t: Bu chhave et al.. 2010). 
Among these, HAT-P-16b has a shorter period, also 
an eccentric orbit, and a much larger radius (1.29 .Rj). 
HD 80606b, on the other hand, has a similar radius, 
and orbits on an extremely eccentric (e — 0.93) orbit 
at 111 day period. It appears that HAT-P-21b is thus an 
unusual, short period, eccentric, massive and compact 
planet. 

The only models from iFortnev et al.l (|2007[) consistent 
with the observed radius are 4Gyr models with 100 
core mass, yielding 1.05 .Rj radius. Probably HAT-P-21b 
has an even higher metal content. HAT-P-21b has a very 
high mean density; it is 8th among all TEPs, and very 
similar to HD 80606b and WASP-14b. 



The flux received by the planet varies between (10.0 ± 
1.5) -10^ ergs^icm~2 and (3.96±0.66) • 10^ ergs'^ cm'^. 
Interestingly, this puts HAT-P-21b on the bordlerline 
between pL (low irradiation) and pM (high irradiation) 
planets. At the time of occultation, HAT-P-21b is just 
approaching its periastron, thus entering the irradiation 
level quoted for pM type planets. 

4.3. HAT-P-22h 

HAT-P-22b has a mass of Mp = 2.147 ± 0.061 Mj, 
radius of Rp = 1.080 ± 0.058 .Rj, and mean density of 
Pp = 2.1ll°-^[|gcm-3. HAT-P-22b orbits a fairly metal 
rich ([Fe/H]' = +0.24 ± 0.08), bright (V=9.732), and 
close-by (82 ± 3 pc) star. Similarly to HAT-P-20, the 
host star has a faint and red neighbor at 9" separation 
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that is co-moving with HAT-P-22 (based on the POSS 
plates and recent Keck/HIRES snapshots), thus they are 
hkely to form a physical pair. 

HAT-P-22b belongs to the moderately massive 
(~ 2 Afj) and compact {Rp w 1 Rj) hot Ju piters, such as 
HAT -P-15b (Mp=1.95Mj, i?p=1.07i?j; iKovacs et afl 
2OT0l). HAT-P-1 4b (Mp=2.23Mj, i?p=1.15i?j; 
Torres et all IMol). and WA SP-8b (Mp=2.25Mj, 
i?p=1.05i?j; iQueloz et aLll2010D . The radius distribu- 
tion is almost bi-modal for these planets (see Figure [T6|). 
with members of the inflated {Rp ~ 1.3 group being: 
HAT-P-23b (Mp=2.09Mj, Rr,=1.37 R.i; this work) 
Kepl er -5b (M„=2 l OMj j?„=1.31 Ry, iKipping fc Bakog 
[2010t IKochetal.1 [2010l). CoRoT-llb (Mp=2.33Mj, 
Rj.=lA3Rr. IGandolfill2010h . 

HAT-P-22b i s bro adly consistent with the models of 
iFortnev et all (poM ). For 300 Myr, 1 Gyr and 4Gyr 
models it requires a 100 M^, 50 and 25 Mq core, 
respectively, to have a radius of ^ 1.08 The low in- 
coming flux (see Table [M]) means that HAT-P-22b is a 
pL class planet. HAT-P-22b can harbor a 0.96 mass 
retrograde moon, which would cause transit timing vari- 
ations (TTVs) of ~ 2 seconds. 

4.4. HAT-P-23b 

HAT-P-23b belongs to the inflated group of 2 Mj plan- 
ets (see discussion above for HAT-P-22b). This planet 
has a mass of Mp = 2.090±0.111 Afj, radius Rp = 1.368± 
0.090 and mean density pp = 1.01 ± 0.18 gcm^^. 
The orbit is nearly circular, with the eccentricity be- 
ing marginally significant. The reason for the somewhat 
higher than usual errors in the RV parameters is the high 
jitter of the star (34.7ms~^), which may be related to the 
moderately high vsini = 8.1 ± 0. 5kms~^ and the very 
close-in orbit of HAT-P-23b. The IFortnev et all (|2008( ) 
models can not reproduce the observed radius of HAT-P- 
23b; even for the youngest, (300 Myr) core-less models, 
the theoretical radius for its mass is 1.25 HAT-P-23b 
orbits its host star on a very close-in orbit. The orbital 
period is only 1.2129 days; almost identical to that of 
OGLE-TR-56b (1.21192days). The nominal planetary 
radius of the two objects is also the same within 1%, 
but OGLE-TR-56b is much less massive (1.39 Mj). The 
flux falling on HAT-P-23b from its host star is one of the 
highest (i.e. belongs to the pM class objects), and is sim- 
ilar to that of HAT-P-7b and OGLE-TR-56b. We also 
calculated th e spiral in-fa ll timcscalc for each new discov - 
ery based onlLevrardrWinisdo erffc r. &: Chabrieil (|2009f ) 
and iDobbs-Dixon et al.l (|2004f ). By assuming that the 
stellar dissipation factor is = 10^, the infall time for 
HAT-P-23b is Onfall = T-^tfl Myr, one of the shortest 
among exoplanets. 

The Rossiter-McLaughlin effect for HAT-P-23b should 
be quite significant, given the moderately high vsini = 
8.1 ± 0.5kms~^ of the host star, and the Ai — 17mmag 
deep transit. The impact parameter is also "ideal" 
(6 = 0.324l:o °oi)' i-e- it is not equatorial (b = 0), where 
there would be a strong degeneracy between the stel- 
lar rotational velocity v sin i and the sky-projected angle 
of the stellar spin axis and the orbital normal. A, and 



is also far from grazing, where the transit is short, and 
other system parameters have lower accuracy. The effec- 
tive temperature of the star {Tcs^, = 5905 ±80 if) is close 
to the critical tem perature of 6250 K noted recently by 
|W mn "eTall (poTol ). which may be a border-line between 
systems where the stellar spin axes and planetary orbital 
normals are preferentially aligned (Teff* < 6250 ii') and 
those that are misaligned (Teg* > 6250 iC) . An alter- 
native hypothesis has been brought up by iSchlaufmanI 
(|2010( ). where misaligned stellar spin axes and orbital 
normals are related to the mass of the host star. The 
mass of HAT-P-23 (1.13 ± 0.04) is sufficiently close to 
the suggested dividing line of Mi, = 1.2 Mq, thus it will 
provide an excellent additional test for these ideas. 

4.5. Summary 

We presented the discovery of four new massive tran- 
siting planets, and provided accurate characterization 
of the host star and planetary parameters. These 4 
new systems are very diverse, and significantly expand 
the sample of ~ 13 other massive (Mp > 2 Afj) plan- 
ets. Two of the new discoveries orbit stars that have 
fainter, most probably physically associated companions. 
The new discoveries do not tend to enhance the mass- 
eccentricity correlation, since only one (HAT-P-21b) is 
significantly eccentric. Also, the tentative mass-wsini 
correlation noted in the Introduction is weakened by the 
new discoveries. The heavier mass planets (HAT-P-20b 
and HAT-P-21b) seem to be inconsistent with current 
theoretical models in that they are too dense, and would 
require a huge core (or metal content) to have such small 
radii. One planet (HAT-P-23b) is also inconsistent with 
the models (unless we assume that the planet is very 
young), but in the opposite s ense; it has an in fiated 
radius. It h as be en noted by I Winn et"al] (|2010D and 
ISchlaufmanI (|2010D that systems exhibiting stellar spin 
axis-planetary orbital normal misalignment are prefer- 
entially eccentric and heavy mass planets (in addition 
to the key parameter being the effective temperature or 
mass of the host star, respectively). The four new planets 
presented in this work will provide additional important 
tests for checking these conjectures. The host stars are 
aU bright (9.7 < V < 12.4), and thus enable in-depth 
future characterization of these systems. 
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